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Abstract 

We study lepton flavor violating (LFV) r and B decays in models with heavy neutrinos to 
constrain the mixing matrix parameters UrN- We find that the best current constraints when 
the heavy neutrinos are purely left-handed come from LFV radiative r decay modes. To obtain 
competitive constraints in LFV B decay it is necessary to probe b XgT^e^ at the 10~^ level. 
When the heavy neutrinos have both left and right-handed couplings, the mixing parameters can 
be constrained by studying LFV B decay modes and LFV r decay into three charged leptons. We 
find that the branching ratios B{t^ ^1^2^3)1 ^{Bs ^ r^e^) and B{h ^ Xsifi^) need to be 
probed at the lO"*^ level in order to constrain the mixing parameters beyond what is known from 
unitarity. 
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I. INTRODUCTION 



In this paper we study lepton flavor violating r and B decays. In the minimal Standard 
Model(SM), generation lepton number is conserved. However, the observation of neutrino 
oscillations implies that family lepton number must be violated [Ij. At present it is not 
clear if the total lepton number is violated. The neutrino oscillation is due to a mismatch 
between the weak and mass eigenstates of neutrinos. This mismatch causes mixing between 
different generations of leptons in the charged current interaction with the W boson. In 
principle, flavor changing neutral current (FCNC) processes in the lepton sector occur as 
well. Some examples would be r — £7, r ^ ^1^2^31 B iF and B il'Xg. Although 
no direct experimental evidence for such FCNC exists, there are experimental constraints 



no direct experimental evic 



12l |. The decays r yu(e)7 have recently been the subject 
of considerable attention [ij"^ They have been studied in connection with LEV occurring 
through mixing with heavy neutrinos in the context of supersymmetric theories where these 
modes are found to be a promising tool to constrain the models. In this paper we investigate 
the potential rates for these processes in Left-Right (LR) models with heavy neutrinos. 

FCNC in the lepton sector that are solely due to mixing in the charged current interaction 
with the usual left-handed W boson and light neutrinos are extremely small because they 
are suppressed by powers of m^/M^ jl^ . One way to increase the FCNC interaction in 
the lepton sector is to introduce heavy neutrinos so that the suppression factor m^/M^ is 
not in effect. This can be done, for example, by introducing a heavy fourth generation. 
If one insists on having just three light left-handed neutrinos, one needs to give the right- 
handed neutrinos heavy Majorana masses. The heavy neutrino can appear through mixing 
in the charged current interaction and enhance the FCNC interaction in the lepton sector. 
The introduction of right-handed neutrinos also raises the possibility of having right-handed 
charged currents by adding to the theory a right-handed W boson. This new charged 
current interaction can lead to additional effects in the above decay processes and here we 
consider such a possibility [l^. 

A natural model of this type is the LR model based on the SU{3)c x SU{2)l x SU (2)^ x 
U{1)b-l gauge group 3]. In the Left-Right model, FCNC interactions arise from several 
sources. In this paper we consider the exchange of Wi^pt bosons at one loop level. It is well 
known that to obtain gauge invariant results one must also include charged Higgs boson 
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effects at one loop level as well as tree level exchange of neutral Higgs bosons We 
will comment on these effects but will not discuss them in detail as they depend on several 
unknown parameters. We concentrate on the effects that depend only on the W mass and 
ignore those that depend on Higgs boson masses to illustrate the constraints that can be 
placed on the mixing with heavy Majorana neutrinos by the processes r — /7, r — > £±£2^3, 
B ££' ^ or 5 — i> ££'Xs. Our paper complements existing studies of the modes /i ^ 67 and 
Kl yue and extends them to include the mixing parameters UrN- 

In LR models there are left-handed light neutrinos and right-handed heavy neutrinos 
i/^, and these neutrinos can be Majorana particles. If there are three light and N heavy 
neutrinos, the general mass term for the neutrinos can be written as 



I 



V 



L 



+ H.C. (1) 
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M'^ is a symmetric matrix which can be diagonalized 

U^M'^U = M'^, (2) 

with the aid of a unitary matrix fj resulting in = diag{mi, m2, m^, M4, M5, ■ ■ ■) with rrii 
and Mi denoting the light and heavy mass eigenvalues respectively. 

If there is no right-handed W-boson interaction, the number of right-handed heavy neutri- 
nos is unrelated to the number of charged leptons. However, when a right-handed W-boson 
is introduced and the heavy neutrinos are required to interact with it, it is natural to have 
the heavy neutrinos and the right-handed charged leptons form SU{2)ji doublets. In this 
case there are as many heavy neutrinos as charged leptons (three). 

The most general charged current interactions of charged leptons and neutrinos with 
W-bosons can be parameterized in the weak interaction basis, as 

9l 



Ciepton = -^W^Lj,{giPL + g'nPR 

Cquark = -^W^U^,[giPL + giPn) D 

_ ^W'^U^, {giPL + -giPa) D, (3) 

where L = (e, fi, r)^, u = [ue, z/^, UrV, z/' = {u'^, z/^, u'^)'^, U = [u, c, t)'^, and 
D = {d, s, 6)^. In the above W and W denote the mass eigenstates of W-bosons with W 
being mostly left-handed and W being mostly right-handed. 



In the mass eigenstate basis, we have [1 



lepton 



^W^Y: ^7. {glPLU,"^* + 9'rPrU, 



R\ m 



1=1 



^W'^U'^j, {giPLV^ + -giPRVn) D^. (4) 

The left-handed and right-handed charged leptons are diagonahzed by the matrices S^'^: 
= S^Il and £^ = S^£r and we have defined the matrices U^* = ELi Sjf'Utj and 
= I]f'=i Sl^'U(^i'+3)j with i = e, fi and r . 

In what follows we will drop the superscript "m" from the fermion fields and always refer 
to mass eigenstates. Note that U^'^ are 3x6 matrices and we construct the matrix 

/ 



tR 




with 



u 



R 



u' 



TI^ TI^ TI^ \ 

'^el ^e2 '^eS '^e4 '^e5 "^66 

TJ^ TJ^ TJ^ JJ^ TJ^ JJ^ 

^ ^ll '^ij,2 '^fii "^m6 

\ I 

^el f/J U^^ 

U^i u% u% u^, 



(5) 



and 



(6) 



I TjR TjR TjR TjR TjR TjR 
\ ^t1 '^t2 ^t3 '^t4 '^tS '^t6 / 

This can be viewed as the unitary matrix which diagonalizes the neutrino mass matrix in 
the basis where the charged lepton mass matrix is already diagonal. The following relations 
hold 



y: u^up^ = 0, 



E ^.?^.'+ E UtuS = S^r (7) 

'■=e,fi,T £=e,fi,T 

There is some information on the matrix elements f/^ and from neutrino oscillation 



experiments Q], which prefer them to be in the ranges 0.50 ~ 0.69 and 0.60 ~ 0.80 respec- 
tively. For the processes that we discuss in this paper there can only be large effects if there 



is substantial mixing with the heavy neutrinos and there are few constraints on these pa- 
rameters. In see-saw models, the generic size of the matrix elements ?7^^(4 5 g) and ^/^^(i 2,3) is 
of order m£,/Mj\f where m^) is the Dirac neutrino mass and is the heavy Majorana neu- 
trino mass. In such a scenario the light neutrino masses are typically m|,/Mjv. This requires 
the heavy neutrino masses to be heavier than a few hundred GeV, and results in the above 
mixing matrix elements being extremely small. In these models, contributions from one W 
exchange to radiatively induced penguin processes, and from box diagrams with two W^s or 
one W and one W exchanges are too small to be observed. Even in these models, however, 
the matrix elements U^^^ 5 can be of order one. For this type of model the exchange of 
W'^s in both radiative penguin and box processes may produce observable effects. There are 
also special cases in which the elements 5 g-^ and 2 3) can be sizeable. An example 
has been discussed in Ref. Iisj, where some of these matrix elements are of order mo/MM, 
but the light neutrino masses, at tree level, are not directly related to them and the ratio 
rrio/MN does not need to be very small. Since our aim is not to study specific models, but 
to provide an estimate of the sensitivity needed in LFV r and B decay modes in order to 
constrain the general Left-Right model mixing beyond the requirement of unitarity, we will 
treat the matrix elements in U' as arbitrary in this paper. 

The couplings g^'^p, and ^^'"^ are in general complex numbers. In renormalizable models 
without left-right gauge boson mixing, g^£'^ = 1, g^jf = 0, g^j]'^ = and ^^'^ = 1. If there 
is left-right gauge boson mixing with a mixing angle {W = Wl cos + Wrsiii^w and 
W = - Wl singly + Wr cosC,w) then 

9l = 9l = cosily, 9^ = 9r = — sin^w] 

9l 

9l = 9l = -— singly, gR = gR = cos^w (8) 
9r 

Throughout the paper we will use the notation: 

J-y^w ^^^w 9l 9l 

For our loop calculations we will assume that all fermions are massless except for the 
top-quark and the heavy right-handed neutrinos. We will also assume that j3 is smaller 
than a few percent in keeping with bounds on W bosons ^|. Finally, we will assume that 
Vt — Wr mixing is small as indicated hj b ^ 57. In particular, following we found in 
20| at the 2a level that there are two allowed ranges for ^g. They correspond to destructive 



and constructive interference with the standard model amphtude respectively and are 

- 0.032 < < -0.027, 

-0.0016 < &g < 0.0037. (10) 

tb 

Overall, C,g is constrained to be very small and in our analysis we will only keep terms linear 
in C,g in the matrix elements. 

We will use this framework to examine LFV induced by neutrino mixing (with new heavy 
neutrinos). Our purpose is to provide an estimate of the sensitivity needed in LFV r and B 
decay modes in order to constrain this scenario beyond the requirement of unitarity of the 
matrix U'. 



II. LFV RADIATIVE r DECAY 

Beginning with jj, ^ e'j j2l|, processes of the form i' £7 have been used to constrain 
new physics, including heavy neutrinos. Here we consider the case of r decay. The one-loop 
effective operator can be calculated in unitary gauge from the diagrams in Figure H We 




w 




I' N I I' N I 



FIG. 1: Diagrams giving rise to i' £7 in unitary gauge. 

first consider the case of a very heavy W so that only the W is exchanged in the loop. 
In keeping with current experimental constraints Eq. we work only to first order in the 
Wl — Wr mixing parameter C,w- This results in two operators which we write in the form 

6 



9l ^ ' 



(11) 



Here F^" is the electromagnetic field strength tensor and the Inami-Lim functions are given 
by 



3A^logA7v ^ 2A^ + - Xn 



4(1-A^)4 

SA^logA^ Ai-llAjv + 4 



;i-A^)3 



2(1-A^)3 4(1 -A^)^ 



(12) 



The exchange of a W can be easily included and it leads to similar expressions: 

, Gf e 

v/2l67r2^ ^ 



C = 4/3- 



9r ^ ' 



(13) 



We now calculate the branching ratio for £' — > £7 neglecting m^. We first consider the 
case without Wj, — Wb. mixing, dominated by the first operator in Eq. El We find 



r(£' ^ £7) 



327r^ 



-VTLn, 



N, 



N 



For r decay, it is convenient to express this result as a fraction of the rate 
to obtain for i = fi, e. 



1927r3 ' 



Rf 



T(t- 



h) 



Tir- 



6a 



N 



(14) 



(15) 



(16) 



The current experimental bounds B(t fi'j) < 3.1 x 10 ^ j^, B(t 67) < 2.7 x 10 ^ 
0], imply that 



Rf, < 1.8 X 10-\ i?e < 1.6 X 10^ 



(17) 



and these in turn, can be used to place the constraints 



N, 



N 



< 1.2 X lo~^ 



< 1.0 X 10" 



:i8) 



If the mixing angles are such that only one heavy neutrino is important, we may use F{x) 
— l/4asx— >ooto estimate that 



ttL* tjL 



< 0.044, 



jtL*tjL 



< 0.13. 



(19) 



For comparison, the experimental bound B{fi — > 67) < 1.2 x 10 jl7[ leads to the constraint 



TjL*TjL 



Similarly, for W exchange one obtains 



< 1.2 X 10"^ 



(20) 



r(f ^ £7) 



327r 



N 



(21) 



Considering once again the case where only one heavy neutrino comes into play, and with 
Aat — > 00 this leads to 



TJ-R TjR* 



< 



0.044 



jjR jjR* 



< 



0.13 

T' 



(22) 



For a typical (3 ~ 0.01 these limits are about an order of magnitude worse than the unitarity 
constraints. The limits become weaker by a factor of four for ~ M^. 

If the mixing parameter is not zero, the second operator in Eq. ^2 can dominate the 
rate because it is not proportional to the light lepton mass. In this case we obtain 



T{i' £7) 



3271^ '' "^^^ 



N 



+ 



N 



N, 



(23) 



With only one heavy neutrino, and using F{x) ^ — 1/4 as a; — »• 00, the constraints are 



-2 ■ 



ttR 



u. 



R 

eN 



U. 



N 

L 

■N 



+ 



+ 



u. 



fiN 



eN 



N 



u. 



eN 



TjR 



< 9.3 X 10" 



< 8.0 X 10" 



< 2.4 X 10^ 



-14 



(24) 



The last result follows from the corresponding analysis for fi —>■ ej. For LR models with 
Wl — Wpt mixing, LFV B decay modes are proportional to C,g making Eq. |^ the most 
stringent constraint in this case. 

III. LFV r DECAY INTO THREE CHARGED LEPTONS 



The pure radiative decays discussed so far do not constrain LR models without Wl — Wr 
mixing because there is no W^W^'y vertex. Similarly, there is no W^W^Z vertex, and the 



modes r 



'1 ^2 ^3 



proceed through box diagrams at leading order. We now derive the 



constraints that can be placed on the neutrino mixing matrix from these modes. 

The effective operator responsible for these decay modes can be calculated from the 
diagram in Figure |21 plus two other diagrams obtained by interchanging Wl ^ Wr and by 
interchanging £i ^ £2- Using dimensional regularization we find 



N 



wvw 

w 



R 



FIG. 2: Box diagram responsible for the decays r ^ £^£2 ™ models with no Wl — Wji mixing. 



a 



27r sin^ 9 



+ u^;j.pt%u^.NUi^,hi.i,PLkhi,i.^^^^ + (^1-^2)}, 



(25) 



S T 

where the Inami-Lim functions E^'j^ are given by 



+ 



+ 1 + 



log/? 



+ 



ElR{XN,,^N,,f3) — \/Xn,Xn. 

+ 
+ 



(A^^,-A;vJ(l-A;vJ(l-/?A^J 

P>^N, log ^jVj 

(Ajv, -A;v,)(l-AArJ(l-/3A^J 
Xn, log Aat^ 



(A^^.-A^J(l-A^J(l-/3A;vJ V 
Ajv, log a 



(A^,-A;v,)(l-A;v,)(l-/5A^,) 
(3/5-1) log/5 



4(l-/5)(l-/?A;vJ(l-/5A^,). 
This result is divergent and we have regulated the divergence by defining 

1 _ 2 

e 4 — n 



+ log47r — 7. 



(26) 



(27) 
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The divergence arises because in LR models these box diagrams are not the only ones that 
contribute to this process. In particular, these models require the existence of neutral scalars 
with tree-level flavor changing couplings . These scalars give rise to a tree- level amplitude 
for this process as well as to several additional one-loop diagrams. To keep our analysis as 
model independent and simple as possible, we will not specify the scalar sector of the left- 
right models. Instead we will use Eq. 1211 drop the 1/e pole, and take a scale /x ~ 1 TeV. This 
approach can be considered as a limit in which the scalars that make the left-right model 
renormalizable are very heavy. To gain some insight into our prescription, we compare our 
result to the complete calculation of Ref. j3| for Kl — > /ie in the appendix. 
Taking both the muon and electron to be massless, we can calculate the rates: 



r r 



'1 ^2 ^3 



1927r3 \^ 1287r2 sin^ 

2 

S 



{ 



NiNj 



+ 64 



LR 



+ 8Re 



^rNi ^hNi Nj ^hNj ^LR 
2 



LR 



EttR* TjL* ttR ttL rpS 
^TNi ^liNi ^l2Ni ^iaN, ^LR 




+ 64 



EttR* ttL* ttR ttL rpT 
^tN^^ liN^^ i2N^^ i^N.-^LR 



EttR* ttL* ttR ttL rpS ttR ttL ttR* ttL* rpT 

^ TNi*^ iiNi'^ l2Nj^ iaNj^LR ^ rNi^ tiNi'^ i2Nj'^ iaNj^ LR 

+ iil^i2)}. 



(28) 



Before comparing with experiment it is convenient to define the ratio 

T(t 



Rl23{T 



'1 ^2 '^3 



'1 ^2 ^3 



(29) 



r(r -> /i UrUf,) 

and normalize the rates this way. When only one heavy neutrino N is important, this 
simplifies to 



L |2 



+ \UrNU^2N\'^\UE*NUf*N\'^ + I f^WV^^^^Af H ^^ff Af f^£^iV H Flr{\n,(3) 



where we have defined the form factor 



Flr{\n, P) 



(30) 



1287r2 sin^ 9w 
+ 64^J^(A^,A^,/3)2 



Flr{^n, Aat, (3) + 8Ej^j^{Xn, Aat, (3)E^^{\n, Aat, (3) 



(31) 
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With the experimental hmits from Ref. 



we find 



R 



liee 



R 



FLRiXN,/3)i 

Flr{\n,P){ 
2Flr{\n,P) 



U. 



tN 



TjR 



U: 



R 

R, 



2Flr{Xn.I3) 



N 



U. 



R 

eN 



+ 



+ 



u. 



R 
tN 



Rp.e.p. 



Flr{Xn, P) 
Flr{Xn, P) 



^ IjlN 



TjR 
TjR 
TjR 

jjR 



U. 



Ut 



N 



u. 



eN 



U. 



u: 



N 



u. 



tN 



U. 



tN 



eN 

TjR 
TjR 

^ fiN 

TjR 

^eN 



)( 

+ 
+ 



'-^fiN 
^^iN 



u 



R 

eN 



eN 



U. 



U. 



u. 



R 

eN 

L 

eN 



+ 



u. 



eN 



TjR 
TjR 



TjR 



N 



+ 



+ 



TjR 

^tN 



TjR 



tN 



'-^/j.N 
^ ^iN 
^eN 



< 1.2 X 10" 



< 1.2 X 10" 



< 1.2 X 10" 



< 2.0 X 10" 



< 1.2 X 10"^ 

< 1.1 X 10"^ 



(32) 



In Figure ITTTl we show the value of FLji{X^,P) for = a.s a function of (3 for selected 
values of Xn- This figure indicates that for a wide range of parameters, -^^^(AAr, /3) is between 



10 



g~n I I r 



"1 I I r~i 




1^ = 100 
^N = 40 
^N = 20 



0.02 

P 



FIG. 3: Flr{Xn,(3) with qr = gi as a function of (3 and \n- 



10^'' and 10^^. With this range in Eq. 1221 we can see that the present experimental limits 
on i?i23 do not yield constraints on the mixing parameters that are significantly better than 
the unitarity bounds. 



A similar exercise for muon decay yields 



eee 



F(/i 



^eN 



TjR 



'-^fiN 



TjR 



+ 



TjR 

'-'/iN 



< 1.0 X 10^1^(33) 
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In this case, if we assume that all the angles are of the same order, we obtain the constraint 
Ulj^ < 0.2. The T decay modes are still far from achieving this level of sensitivity. If we 
write for example. 



R 



liee 



TjL* 



jjL* 



eN 



f/5 



tN 



eN 



r(/i 



eee 



r(/i 



< 1.0 X IQ-^^ 



tN 



TjL* 



+ 



R 
tN 



(34) 



we see that the r decay bounds need to improve by six orders of magnitude in order to 
be competitive with the already available muon decay limit. Of course, the muon decay is 
not sensitive to all the parameters needed to describe neutrino mixing in general and the r 
decay data is complementary. 

IV. LFV B DECAY: OPERATORS 



We now turn our attention to B decay modes and start by calculating the basic quark 
level LFV process h djii'. In unitary gauge the process occurs through the box diagram 
of Figure m. We distinguish several CcLSGS cLS before: left-handed heavy neutrinos; LR models 



b 



w 

i/W\A( 



>k N 



l/WW 

w 



FIG. 4: Box diagram responsible for the process b djU'. 

with Wl — Wr mixing; right-handed heavy neutrinos; and LR models without Wl — Wr 
mixing. 

A. Left-handed heavy neutrinos 



A straightforward calculation of the diagram in Figure ID produces the operator: 



72 Stt sin^ e 



w 



(35) 
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where dj refers to a c? or an s quark. With the aid of the unitarity relations of Eq. [71 we find 

3 



N 



(1-A^)(l-Ai) 



{4-8Xt + Xl)logXt (4-8A;v + A^)logAjv | 
(A^-Ai)(l-At)2 + (A,-A;v)(l-A^)M ■ ^ ^ 

This result is in agreement with the existing result for Kl — > /x^e^ [l5| when the b quark is 

replaced by an s quark. 



B. Left-right model with Wl — Wr mixing 

In LR models with mixing there is a second operator that can be obtained from Figure U 
^27rsin2 Ow ^ 

+ EliK, A;v)^7.7m [u^nU^^Pr + U^nU^nPl) i'd.YY {Vt^^ViPn + V,^;V,^,Pl) b . 

(37) 

The Inami-Lim functions are the same as those in Eq. |2niwith (5 = 1. As mentioned above, 
these operators produce observables proportional to and cannot place constraints that 
are competitive with LEV radiative r decay. 



C. Right-handed heavy neutrino 

Models with a mostly right handed heavy neutrino would proceed through the diagram 
in Eigure m with two W bosons exchanged. Ignoring the Wl — Wr mixing this results in an 
operator 

^ = "2. V,^;V,^U^^U^*Mf3Xu (3XN)h,PRi'd,YPRb (38) 

A glance at Eq. EEl reveals that this operator will make contributions to LEV B decay rates 
that are suppressed by at least a factor of (3^. 



D. Left-right models without Wl — Wr mixing 



In this scenario the heavy neutrinos have left and right-handed couplings and both the 
W and W' appear. The LEV operator for B decay arises from box diagrams like the one in 

13 



Figure ^ with one W and one W, we find 

^ = ^ 27rsin^g^ ^^ [El^{Xt,XN,P) {v,'^:Y,^U§,U^,^lPd' d.Pnh 

+ l^,^f/ii^f/,?^^7.7M^H^'rf,7M7.i^L&)l , (39) 



where the Inami-Lim functions were given in Eq. 

For the process B ££' the two operators in Eq. EHlcan be reduced to one by using the 
relation 

'lu7t^®'^y^^7u 4(1 O 1) + cr^i, (g) cr^^, (40) 
and dropping the last term in anticipation of the vanishing of the matrix element 

<0|J,-a^,(l±75)6|5j>=0. (41) 

We obtain, 



'■5 



y Z ^TT Sm U]y ^ ' 

(y,^a*ViU^j,Ui^,^IPd'diPRh + V,^;V,lU^;U^*^IPd'd,PLh) . (42) 



X 



V. LFV B DECAY PHENOMENOLOGY 



We now use the operators obtained in the previous section to compute their contribution 
to selected LFV B decay modes. 



dd 



We first consider the mode B r^£^, the analogue of the Ki^ — > [de^ mode which has 



been discussed extensively in the literature 



15 



For a heavy left-handed neutrino we find after summing the two modes and neglecting 
the mass of £ = /i, e, 



r(5, ^ dd) 



a 



256 vr Uvrsin^^ 



w 



FlmlMu 1 



\yL*yL\2 



N 



(43) 
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When the heavy neutrino has both left and right handed couphngs (with an exphcit W 
and no mixing) we find 



32 71 V47rsin2^ 



w 



^"-"B n2 
B 2 /^O 



/5 1 



+ 



AT 

E • f/i^f/.'^; \ElniK Aiv, /3) + 4EL(A*, A^, /?) 



AT 



(44) 



For numerical purposes it is natural to compare these rates to the standard model rate 
for Bd, -^T+T~, 



r(5, 



T T 



a 



TT V 47r sin^ 6'vi/ 

where the Inami-Lim function Y{Xt) ~ 1.06 [2^. We define 

r(5, - r±f^) 



\v:,vt,\'Y\\t 



(45) 



r(5,^r+r-)' 

and in terms of this definition we find for left handed heavy neutrinos 



(46) 



RBri = 3.7 X 10" 



AT 



(47) 



To compare the sensitivity of the modes B t^£^ and r ^ £7 to the neutrino mixing 
parameters we plot in Figure El the ratio {El/ FY as a function of Xn- 

Assuming there is only one heavy neutrino and using Eq. UHlas well as {El/ F)"^ ~ 7500 
from Figure El we find 



Rbt), < 3.3 X 10~^ RBre < 2.8 X 10"^ 



(48) 



The standard model expectations for Bd t^t j23| with the central values for the CKM 



angles found in 



2J| are 



B{Bs T+T-) = 1.1 X 10"^ B{Bd T+T-) = 3.3 x 10"^ 



(49) 



Consequently, one would need a single event sensitivity of at least 10~^ for B{Bs — > r^e^) 
(10^^ for B{Bs r^yU^)) to improve on the existing constraints from radiative r decay. A 
glance at Table U indicates that one would need an order of magnitude improvement over 



15 



10000 



UJ 



7500 
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20 



40 



FIG. 5: Ratio of form factors {Ei/F)"^ as a function of A 



AC- 



TABLE I: Summary of current experimental bounds for B l^£' 





Branching Ratio 


B e±/uT 
B e^T^ 

B ll^T^ 


< 1.8 X 10-'^ Babar [7j 

< 1.7 X lO-'^ Belle 

< 1.1 X 10'^ Cleo [2] 

< 3.8 X 10-5 Cleo [2] 



the current best Umit from Belle for B — > e^/i"^. There is some hope that this sensitivity 
may be attainable in the future. For example, the estimated single event sensitivity for the 
B ^ HH modes at CDF with 15 fb~^ is j25[ 



1.3 X 10"^ for Bs ^^^^ 
4.7 X 10"^° for Bd ^^n^. 



(50) 



Of course, there are additional experimental difficulties for modes involving a r lepton, but 
we may regard Eq. EHl as the benchmark needed to improve upon limits from radiative r 
decay. For comparison, the same analysis applied to the limit B{Kl — > /x^e^) < 4.7 x 10^^^ 
jl7|. yields the bound 



TjL*TjL 



< 0.07, 



(51) 
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which is also weaker than the corresponding bound from /i — 67, Eq. 1201 
For right handed heavy neutrinos, Eq. 07] becomes 



3.7 X 10 



-3 



N 



(52) 



The bounds obtained are worse than those for left-handed heavy neutrinos. 

Finally for the case where the heavy neutrino has both left and right handed couplings 
and the Wl — Wr mixing can be ignored we obtain 



R 



Brl 



0.36/3; 



9 |T/*T/ 12 



N 



+ 



N 



(53) 



It is harder to interpret this result because there are several unknown parameters. To gain 

1 



I I r 



uj 10 



"I I I r~i 




A,fj = 100 
^N = 40 
^N = 20 



0.02 

(3 



FIG. 6: E'IjI as a function of /? for selected values of Aat. 



some insight into this result we consider the simplified case with only one heavy neutrino 
discussed in the introduction. We further assume that V^^ ~ V^^, and that = Qr- We 
define the form factor 

2 



^ Q.36p-' El^{X,, Xn, (3) + AEij,{Xu A^, (3) 



(54) 
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and plot it in Figure 6 as a function of (3 for selected values of the heavy neutrino mass. We 
see from the figure that Ej^jl can be between 0.01 and 0.1 for a wide range of parameters. 
Combining this result with Eq. [321 we find, 



RBri ~ (0.01 - 0.1) llu.Wj,]'' + \u[;^u^^\^' 



(55) 



This in turn implies that RbtI has to be probed at the 10~^ level to constrain the neutrino 
mixing parameters beyond what is known from unitarity. Considering the SM expectation 
Eq. |3ni this implies a benchmark number for B{B ^ r^£^) at the 10^^ level, probably 
beyond reach for the foreseeable future. The corresponding bound B{Kl f^^^^) < 
4.7 X 10"^^ [Ui] leads in turn to. 



< 0.03. 



(56) 



It is instructive to compare these modes to the r~ —>■ iii2^t modes. To this effect we 
plot in Figure [Tfa) the ratio of form factors E^^^/Flji as a function of Aat for /3 = 0.0065 
(which corresponds to Mw' ~ 1 TeV). Using e)/^jI ~ lO^^F^/j from Figure (a) we write 




10 



J I I I I I L 



20 



40 



0.3 



liJ 



0.2 



0.1 - 



T — I — r 
(b) 




FIG. 7: Squared ratio of form factors for a) Ej^j^/Flji and b) E^^j^/Flji as a function of A at for 
P = 0.0065. 



'tN^TN\ 



_|_ \ TJ^*TJ^ 



N\ 



(57) 



which can be compared directly to Eq. 
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B. Inclusive modes b — > sr^£^ 

We turn our attention to the inclusive process b sr^i"^. This choice is motivated 
by several factors. A semileptonic mode removes the helicity suppression present in the 
B T^l^ modes for the case of left-handed heavy neutrinos. ^ Requiring strangeness in 
the final state removes the CKM suppression factor Vtd/Vcb- Finally, considering an inclusive 
process removes the QCD suppression factor Fb/Mb- 

For left handed heavy neutrinos we find 

„2 \ 



a 



1927r3 \5127r2sin^^ 



w . 



N 



(58) 



where 

J(x) = 1 - Sx^ + - - 24x^ log(x) (59) 

is the usual kinematic factor for a non-zero r mass and I{mr/mb) ~ 0.33 with mb ~ 4.5 GeV. 
Using a b lifetime given by 



r(6) = 5.8 



1927r^ 



rrib 



\Vcb\^ 



(60) 



with I{mc/rnb) ~ 0.5 (for rric ~ 1.4 GeV), this gives 



B{b djT^i^) ^ 2.6 X 10 



-8 





2 







N 



(61) 



To compare with the radiative r decay modes we assume that there is only one heavy 
neutrino, use Eq. [THlas well as \El/F\^ ~ 7500 from Figure El (a) to obtain. 



B{b- 




< 


2.3 X 10"^ 


B(b- 




< 


2.0 X 10"^ 


B{b- 


^ dr^fi^) 


< 


9.1 X 10"^° 


B{b- 


dT^e^) 


< 


7.7 X 10^^ 



(62) 

These results imply that one needs at least a 10^'' sensitivity in B{b sr'^e^) to obtain 
constraints competitive with the radiative r decay for left-handed heavy neutrinos. Table ITTl 



^ In the kaon sector the hehcity suppression may be removed by considering instead K —> True modes. The 
best bound in that case, \U^^Ul^]^\ < 2, arises from B{K~^ 7r+/x+e~) < 2.8 x 10" 
than the unitarity constraint. 



True 



but is worse 
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indicates that current B-factory results are close to this benchmark, although there are 
significant experimental hurdles for detection of r leptons. 



TABLE II: Summary of current experimental bounds for B Xi^i' modes. 





Branchin 


g Ratio 




b- 




< 2.2 X 10-5 


CLEO 


[lOj 


B - 


> TTe^fi^ 


< 1.6 X 10"^ 


CLEO 


M 


B 




< 4.0 X 10"^ 


BaBar 


[12] 


B - 




< 1.6 X 10"^ 


CLEO 




B - 


^ pe^fi^ 


< 3.2 X 10"^ 


CLEO 




B 




< 3.4 X 10"^ 


BaBar 


[12] 


B 




< 6.2 X 10"^ 


CLEO 





As an estimate for the reach of future experiments we start from the Tevatron studies 
indicating that with 2 fb"^, CDF could detect 61 Bd K*^fifi events assuming a branching 
ratio of 1.5 x 10~^ [2^. We can turn this number into an approximate single event sensitivity 
with 15 fb^^ of 3.2 x 10"^ for this mode. This in turn indicates that improved constraints 
are possible, at least from the b sre mode. 

For heavy neutrinos with left and right handed couplings but vanishing Wl — Wr mixing 
we find 

>,2 \ 



Glml 



a 



1927r3 \ l287r2 sin^ 



rrir 
nib 



X 



s 

LR 



N 



+ 64 



T 

LR 



N 



'^LR 



N 



N 



(63) 



If we assume that only one heavy neutrino is important, that ~ gi, and that V^f ~ V^^, 



we can write 

B(b djT^i^) 

where we have defined 



Vcb 



TjR TjL 



+ 



ttL*ttR* 



ELR.i'^t, Aiv, (3) 



(64) 



£;S(A„ A^., /?) ^ 1.04 X io~'p-' {El^f + %El^Ei^ + 64(£;^^; 



n5 \2 



(65) 



20 




FIG. 8: E^ll as a function of /3 for selected values of A at with gn = ql- 

From Figure |H1 we see that E^^]^ ~ 10"^ for a wide range of parameters. This allows us to 
write 



B(h ^ ST^I^) ~ 1 X 10-^ (yMN\ + \UmU?N\ ) ■ (66) 

A sensitivity of at least 10~^ is thus required for the mode h — > sr^£^ to place significant 
constraints on the neutrino mixing parameters. 

To compare these modes to the r~ ^1^-2^3 modes, we plot in Figure [7||^b) the ratio of 
form factors E^)^/ Fip> as a function of Xn for /3 = 0.0065. Using E^^l ~ 0.02^^^ we write 



B{b sT^i^) ~ 0.02Flr 



TjR jjL 



+ 



TjL*TjR* 



(67) 



which can be compared directly to Eq. |221 



VI. SUMMARY AND CONCLUSIONS 

We have studied LFV r and B decay modes within the context of neutrino mixing with 
additional heavy neutrinos. We have considered generic left-right models and distinguished 
three scenarios. In the first scenario we consider, the W has a negligible effect and the lepton 
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interactions are purely left-handed. In this case the best constraints on the parameters 
describing the neutrino mixing arise from radiative r decay. Assuming one very heavy 
neutrino dominates, the current best bound is placed by r — > /i7, 

\U;:*^U,^M < 0.044. (68) 

To obtain a competitive constraint from Bs — > T^fJ'^, this mode would have to be probed at 
the 10^^ level as indicated by Eqs. |3HlaiidE3 Similarly, to obtain a comparable constraint 
in any of the inclusive modes, b — > sr^e^ needs to be probed at the 10"'' level as indicated 
in Eq. El 

The second scenario we considered involved a very heavy W. In this case the effect of the 
right-handed interaction can only be felt at low energies through Wl — Wr mixing. In this 
scenario the best constraints arise from radiative r decay and lead to unobservably small 
rates in LFV B decay modes. 

The third and final scenario we considered is one in which both the W and the W play 
a role in the lepton charged currents, but there is no Wl — Wr mixing. In this more general 
case, the (6x6) mixing matrix in the neutrino sector has many unknown parameters and 
as a practical matter different decay modes will in general probe different combinations of 
these parameters. As a benchmark for the sensitivity needed to probe this scenario we have 
considered three simplified cases with results summarized in Table UTTl The salient features 
are 

• Radiative r decay modes do not probe this scenario in the limit of no Wl — Wr mixing. 

• Unitarity implies that the matrix elements of the neutrino mixing matrix satisfy 
\U(^'^\ < 1. To place significant constraints (better than the unitarity limit), i?i23 
has to be measured with a sensitivity of at least 10~^, two orders of magnitude better 
than current limits. 

• The benchmark for significant constraints from B decay is a sensitivity of 10~^ for 
B{Bs r±£T) and of 10"^ for B{h sr^l^). 

In all cases the sensitivity to the neutrino mixing parameters is much smaller than what 
already exists from the study of ^ ^ ej and Kl — > fie modes. However, LFV r and B 
decay modes offer an opportunity to complement those results by providing constraints on 
the Utn mixing angles. 
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TABLE III: Summary of results for a heavy neutrino with left and right handed couplings. We 
consider three cases with only one non-zero mixing with a heavy right-handed neutrino: a) 7^ 0, 
b) / 0, and c) U^^ / 0. 
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APPENDIX A: COMPARISON WITH THE Mh ^ 00 LIMIT OF REF 



In this appendix we compare the form factor in M — if £2 (for a spinless meson M) 
as obtained from Eq. |2ni and the complete one-loop calculation of Ref. 3]. In a LR 
model with a simple scalar sector is considered and a complete (gauge independent and 
finite) result for the form factor is obtained. In addition to the gauge boson contribution 
that we consider in this paper, the complete result depends on the parameters of the scalar 
sector, namely at least three scalar masses and scalar sector couplings. To compare this to 



our result we take the Higgs masses to infinity in the result of Ref. 



In what follows we use the notation in Appendix-B of We drop a common factor 



1^ 



containing the mixing angles, the coupling constants and the factor vAA/v- The diagrams 
that contribute to the form factor in the limit when the Higgs masses are very heavy are: 



3a+3b - box diagrams with one W and one W (these are the ones we consider in 
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unitary gauge) 

{FF)a, = (3 



4 



where the functions J2 and Fi are defined in 



(Al) 



3c+3d+3e+3f - vertex corrections involving the exchange of the neutral scalars with 
tree-level fiavor changing couplings. 



(FF) 



cdef 



4 



A^logAi 



+ 



A^ log \n 



(l-Ai)(l-/3Ai) (l-A^)(l-/5A;v) 



+ 



(2/jAtAjv-A,-Ajv)log/? 
(l-/5)(l-/?At)(l-/?A;v), 



+ 



(A2) 



l+3g+3h - tree-level exchange of scalars with fiavor changing couplings as well as 
self-energy corrections to this. 

p /log/5 , , ^ ..f 1 ^ 



(FF) 



gh 



+ logA<i-2 



(A3) 



4 VI -/? 

where A,^ = M|/M^ and is the mass of the neutral scalars (assumed degenerate 
in 13]) with tree- level flavor changing couplings. 

• all other diagrams considered in Ref. yield contributions to the form factor that 
vanish as inverse powers of Mh in the limit of inflnite mass for the physical Higgs's 
that occur. 

The final result in the limit Mh ^ 00 is then 



(FF) 



+ 



(A^-Ai)(l-Ai)(l-/3Ai) 
AivlogAjv(/?A^-/3AAr-Ajv + 4) , 3/3 log /3 



(A4) 



(Ai - A^)(l - Xn){1 - f3XN) (1 -/?)(!- /?Ai)(l - f3XN)_ 
This is to be compared to our result for the same form-factor (which appears in Eq. 
as obtained from Eq. 1211 {FF)us ~ (_Elji{Xt, Xn, P) + 4:E]^f^{Xt, Xn, P)) and normalized in 
the same way as the result of Ref. |,L5| : 

'l / ,,2 \ \ ^^^\ fn\2 



(FF) = f 



2:1 I /^M I . I A,logA,(/3A^-/3A,-A, + 4) 
e+°^lM2,j+ + (A^_A,)(1-A,)(1-/?A,) 



+ 



Xn log Xn {PX% - PXn - Ajv + 4) 



3/5 log /3 



(Xt - A^)(l - Aiv)(l - PXn) (1 - /3)(1 - PXt){l - PXn) 



(A5) 
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Comparing these two results after using our prescription shows that we reproduce the com- 
plete calculation in the limit Mh ^ oo up to the numerical factor 3 for the choice /i = Mh- 
This number depends on the renormalization scheme used and is of the same order as the 
logarithm log(/i^/M^) ~ 5 for the scale yU ~ 1 TeV we choose. 
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